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Interactions  Between  Wetlands 

us  Army  Corps  and  Tidal  Inlets 

of  Engineers® 

by  Alejandro  Sanchez 


PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  presents 
numerical  simulations  investigating  how  the  loss  of  wetlands  in  estuaries  modifies  tidal 
processes  in  inlet  navigation  channels.  The  implications  for  wetland  loss  and  construction  of 
wetlands  to  tidal  inlets  and  navigation  channels  are  discussed,  and  hypotheses  are  presented  for 
possible  engineering  approaches  to  mitigate  wetland  loss  and  improve  the  sedimentary  sustain¬ 
ability  of  constructed  wetlands. 

INTRODUCTION:  One  of  the  navigation  missions  of  the  U.S.Army  Corps  of  Engineers 
(USACE)  is  to  maintain  navigable  waterways  at  coastal  inlets  and  entrances.  This  mission 
requires  an  understanding  of  navigation  channel  shoaling;  environmental  relationships  among  the 
inlet,  estuary,  and  adjacent  beaches;  and  how  modifications  to  the  inlet  and  estuarine  system 
(including  placement  of  dredged  material  to  create,  restore,  or  enhance  wetlands)  may  alter 
channel  shoaling  and  adjacent  morphologic  features.  In  fiscal  year  2007,  the  USACE  dredged 
more  than  181  million  cu  yd  (138  million  cu  m)  of  sediment  to  maintain  navigation  through  inlet, 
inland,  and  estuary  channels  (USACE,  National  Data  Center  2008).  The  majority  of  this 
sediment  is  not  contaminated  and  can  be  accessed  to  create,  restore,  and  enhance  wetlands  under 
appropriate  economic,  logistical,  and  regulatory  constraints.  Placing  dredged  sediment  to  sustain 
environmental  resources  such  as  wetlands  directly  supports  USACE  Environmental  Operating 
Principles  to  “design  solutions  that  support  and  reinforce  one  another”  and  “seek  to  mitigate 
cumulative  impacts  to  the  environment”  (USACE  2008). 

The  USACE  has  been  involved  in  the  creation  and  restoration  of  wetlands  for  more  than  30  years 
through  beneficial  use  of  dredged  material  and  as  mitigation  projects  (Woodhouse  et  al.  1974; 
Garbisch  et  al.  1975;  Landin  1982;  Palermo  1992;  Dunne  et  al.  1998;  Yozzo  et  al.  2004).  New 
and  restored  wetlands  are  usually  designed  by  optimizing  their  functionality  (benefit  or 
performance)  as  determined  by  the  hydrology  hydraulics,  sediment  transport,  water  quality,  and 
wildlife  support.  Wetland  establishment  and  restoration  are  long-term  processes  that  depend  on 
sustained  and  conducive  physical,  chemical,  and  geological  conditions.  These  conditions  are 
influenced  by  estuary-wetland  interactions  such  as  water  storage  and  sediment  pathways.  These 
interactions  have  been  little  explored  in  the  design  and  restoration  of  wetlands  and  in  the  design 
and  functioning  of  inlets  and  navigation  channels.  As  discussed  in  this  CHETN,  the  coverage, 
extent,  and  density  of  vegetation  associated  with  estuarine  wetlands  influence  the  long-term 
evolution  of  estuary  morphology  and  tidal  inlets,  navigation  channels,  and  the  wetlands 
themselves.  This  knowledge  will  be  translated  into  predictive  morphological  tools  for  manage¬ 
ment  of  navigation  channels  and  tidal  inlets  as  part  of  the  Estuaries  and  Channels  Work  Unit  of 
the  Coastal  Inlets  Research  Program  (CIRP). 
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WETLANDS  AND  ESTUARINE  DYNAMICS:  A  leading  influence  of  wetlands  on  estuaries  is 
their  contribution  to  determining  the  tidal  prism  (the  volume  of  water  entering  an  inlet  during 
flood  tide  or  exiting  an  inlet  during  ebb  tide).  As  wetlands  are  lost  or  restored,  the  tidal  prism  is 
modified.  Changes  in  tidal  prism  have  several  implications  on  the  estuarine  hydrodynamics  and, 
consequently,  sediment  transport  and  morphology.  Many  inlets  serve  as  navigation  routes  for 
commercial  and  recreational  vessels.  An  increased  tidal  prism  caused  by  wetland  loss  increases 
current  velocity  through  the  tidal  inlet,  thereby  increasing  the  inlet  cross-sectional  area  and 
sometimes  opening  new  pennanent  or  ephemeral  inlets.  Conversely,  the  development  or 
construction  of  wetlands  within  an  estuary  reduces  bay  area  and  the  tidal  prism,  which  will 
change  the  phase  and  magnitude  of  the  ebb  and  flood  tides.  The  non-uniform  change  in  tidal 
prism  within  the  estuary  can  significantly  alter  the  circulation  pattern  and  create  new  channels 
and  fill  in  existing  ones  (FitzGerald  et  al.  2004).  Wetlands  can  be  restored  by  opening 
(breaching)  diked  sections  of  the  shore.  Breached  restoration  sites  have  been  linked  to  an 
increase  in  estuarine  cross-sectional  area  and  erosion  of  the  existing  salt  marsh  (Pethick  2002). 

Tidal  wetlands  also  influence  estuarine  dynamics  by  producing  tidal  asymmetry.  For  marsh 
platforms  that  lie  above  low  tide  elevation,  the  increased  flow  resistance  by  wetland  vegetation 
during  high  tides  (as  compared  to  no  resistance  at  low  water)  produces  longer,  but  weaker  flood 
tides  and  shorter,  but  stronger  ebb  tide — called  ebb  dominance  (Boon  and  Byme  1981;  Speer  and 
Aubrey  1985;  Friedrichs  and  Aubrey  1988).  By  continuity  of  mass,  the  current  velocity  at  an 
inlet  is  directly  proportional  to  the  rate  at  which  the  tidal  prism  changes  with  water  elevation. 
The  convex  bathymetric  profile  associated  with  wetlands  produces  an  increased  rate  of  change  in 
tidal  prism  during  ebb  tide,  enhancing  ebb  dominance  (Boon  1975;  Dronkers  1986).  Because 
sediment  transport  is  proportional  to  the  current  velocity  to  a  power  greater  than  one,  tidal 
asymmetry  determines  the  net  direction  of  the  depth-averaged  sediment  transport  flux.  Therefore, 
flood-dominant  systems  will  have  a  net  import  of  sediments,  whereas  ebb-dominant  systems  will 
have  a  net  export  of  sediments.  The  influence  of  wetlands  on  tidal  asymmetry  depends  on  the 
extent  of  wetland  coverage,  the  tidal  amplitude,  the  size  of  estuary,  and  other  factors.  Earlier 
studies  of  tidal  asymmetry  in  estuaries  have  consisted  of  analytical  solutions  or  numerical 
simulations  of  simplified  geometries  and  semi-diurnal  tides  (Speer  and  Aubrey  1985;  Friedrichs 
and  Aubrey  1988;  Lawrence  et  al.  2004;  Fortunato  and  Oliveira  2005).  To  improve 
understanding  of  the  interactions  between  wetlands  and  their  inlets  and  channels,  it  is  necessary 
to  study  estuaries  under  mixed  tide  and  with  realistic  bathymetry  and  geometry.  This  CHETN  is 
intended  to  be  a  step  beyond  the  previous  work  in  that  it  applies  a  two-dimensional  (depth- 
averaged)  circulation  model  to  an  idealized  estuary  with  a  bathymetry  and  geometry  that  more 
closely  resemble  the  natural  system. 

During  stonns,  coastal  wetlands  dissipate  wave  energy,  reducing  shore  and  bay  bottom  erosion, 
thus  preserving  the  tidal  prism  and  augmenting  the  functioning  of  shore-protection  structures 
during  storms  (Owen  1984;  Moeller  et  al.  1996).  Wetlands  and  complex  channel  networks  absorb 
tidal  energy  and  reduce  tidal  amplitudes  and  currents.  In  addition,  vegetation  roots  shelter  and 
reinforce  bottom  sediments  to  protect  the  area  from  erosion  (Coops  et  al.  1996).  Wave  setup 
contributes  to  elevated  water  level  during  storms  and  can  account  for  a  considerable  part  of  total 
storm  water  level.  Analytical  work  of  Dean  and  Bender  (2006)  recently  showed  that  vegetation 
can  significantly  reduce  wave  setup  and  even  produce  a  setdown. 
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Wetlands  are  a  significant  part  of  the  regional  sediment  budget,  in  which  they  serve  as  sediment 
sinks  even  in  ebb-dominated  tidal  regimes.  For  example,  Ganju  et  al.  (2005)  estimated  a  net 
import  of  14  metric  tons  over  a  34-day  period  in  two  ebb-dominated  tidal  channels  in  a  small 
marsh  island  in  San  Francisco  Bay,  CA.  Wolaver  et  al.  (1988)  measured  suspended  sediment  flux 
of  827  g/m  / year  into  a  marsh  in  North  Inlet,  SC.  The  influence  of  wetland  loss  on  the  regional 
sediment  distribution  is  poorly  understood.  A  large  portion  of  the  sediment  accretion  in  marshes 
is  produced  in  situ  by  low-ground  production  (mostly  from  vegetation).  In  some  estuaries,  the 
resuspension  and  advection  of  marsh  sediments  are  related  to  episodic  meteorological 
phenomena  such  as  tropical  stonns  and  frontal  passages.  Baumann  et  al.  (1984)  found  that 
70  to  80  percent  of  the  accreted  sediment  in  Barataria  Bay,  LA,  salt  marshes  was  deposited 
during  winter  frontal  passages.  Stumpf  (1983)  also  found  similar  storm-dominated  sediment 
deposition  in  Delaware  Bay  marshes.  It  is  generally  accepted  that,  when  wetlands  erode,  most  of 
the  fine  sediments  are  transported  offshore  and  removed  from  the  littoral  system.  Determining 
sediment  loss  and  gain  in  estuaries  during  typical  and  severe  storms  is  necessary  to  estimate 
long-term  estuarine  evolution. 

Wetland  processes  and  their  interactions  with  tidal  inlets  and  navigation  channels  must  be 
integrated  to  predict  long-term  morphological  evolution.  This  CHETN  illustrates  how  idealized 
numerical  modeling  can  be  used  to  examine  and  quantify  these  interactions  and  reveal  areas 
lacking  knowledge  for  further  model  development. 

LONG-TERM  EXPERIENCE:  Data  spanning  more  than  100  years  in  Barataria  Bay,  LA, 
illustrate  how  wetland  loss  can  modify  tidal  prism,  inlet  area,  and  ultimately  the  regional 
sediment  budget.  Because  of  rapid  relative  sea  level  rise  (1.03  cm/year),  wetland  loss,  and 
erosional  processes  within  Barataria  Bay,  the  bay  area  increased  more  than  775  km  from  1956  to 
1990.  The  resulting  increase  in  bay  area  increased  the  tidal  prism  through  the  four  inlets  fronting 
the  estuary,  and  increased  the  inlet  cross-sectional  areas  and  shoal  volumes.  For  example,  the 
ebb-tidal  shoal  at  Barataria  Pass  moved  more  than  2.2  km  seaward  from  1880  to  1980,  the  inlet 
cross-sectional  area  increased  from  4,200  to  6,900  m  ,  and  the  pass  deepened  from  14  to  37  m. 
Sediment  to  build  the  expanding  ebb  tidal  shoals  originated  from  deepening  of  the  inlets,  the 
longshore  sediment  transport  system,  and  the  adjacent  barrier  islands.  Thus,  wetland  loss  and 
relative  sea  level  rise  ultimately  resulted  in  a  sink  to  the  regional  sediment  budget  (FitzGerald 
etal.  2004). 

In  multiple  inlet-bay  systems,  wetland  loss  can  switch  the  dominant  inlet  and  change  the  regional 
sediment  transport  pattern.  For  example,  the  greater  loss  of  wetlands  in  the  northern  part  of 
Barataria  Bay  directed  a  larger  proportion  of  the  tidal  prism  through  Barataria  Pass  with  respect 
to  Caminada  Pass,  Pass  Abel,  and  Quatre  Bayou  Pass,  and  eventually  realigned  the  channel  to  a 
more  north-south  alignment.  The  change  in  channel  orientation  caused  the  ebb  tidal  shoal  to  shift 
further  east,  eventually  providing  more  protection  to  the  western  Grand  Isle  barrier  island 
(FitzGerald  et  al.  2004).  This  long-tenn  data  set  at  Barataria  Bay  lends  insight  into  the  possible 
role  wetlands  can  play  in  changing  tidal  prism,  inlet  cross-sectional  area,  and  flow  through  inlet 
channels  for  the  benefit  of  navigation  and,  in  this  case,  for  reducing  operation  and  maintenance 
costs.  Due  to  relative  sea  level  rise,  wetland  loss  is  expected  to  increase  in  Barataria  Bay  as  in 
many  other  locations.  The  Barataria-Terrebone  Ecological  Landscape  Spatial  Simulation  model 
predicts  an  almost  90  percent  wetland  loss  in  the  Barataria  Basin  from  1988  to  2018  for  the 
present  rate  of  sea  level  rise  (Martin  et  al.  2002). 
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NUMERICAL  EXPERIMENT:  The  change  in  current  velocity  at  tidal  inlets  caused  by  wetland 
loss  is  examined  here  through  numerical  modeling  with  five  idealized  estuaries  representing 
varying  degrees  of  wetland  coverage.  The  hydrodynamics  are  simulated  with  the  Coastal 
Modeling  System  (CMS)  circulation  model  CMS-Flow  (Buttolph  et  al.  2006)  that  solves  the 
conservative  form  of  the  non-linear  shallow- water  equations.  The  CMS  is  being  developed  under 
the  CIRP  and  can  interactively  simulate  waves,  circulation,  sediment  transport,  and  morphology 
change. 

Idealized  Estuaries.  The  geomorphic  classification  of  estuaries  consists  of  four  main  groups: 
(1)  coastal  plain  or  drowned  river  valley  estuaries  (e.g.,  Chesapeake  Bay,  MD),  (2)  lagoonal  or 
bar-built  estuaries  (e.g.,  Laguna  Madre,  TX),  (3)  fjord-type  (e.g.,  Penobscot  Bay,  ME),  and 
(4)  tectonically  created  estuaries  (e.g.,  San  Francisco  Bay,  CA)  (Pritchard  1967).  This  CHETN 
discusses  lagoonal  or  bar-built  estuaries  (Figure  1).  ft  is  hypothesized  that  these  estuaries  are 
likely  to  be  most  affected  by  relative  sea  level  rise  and  wetland  loss  because  of  their  shallow 
water,  narrow  coastal  inlets,  and  extensive  wetlands. 

The  following  parameters  were  considered  in  the  design  of  the  idealized  estuaries: 
(1)  geomorphic  classification,  (2)  basic  dimensions  of  width,  length,  surface  area  and  depth, 
(3)  ratio  of  tidal  amplitude  to  average  bay  depth  (relative  tidal  amplitude),  (4)  variation  of  the 
water  surface  area  with  depth  (called  the  hypsometry),  (5)  ratio  of  volume  of  water  over  wetland 
to  bay  water  volume,  and  (6)  ratio  of  wetland  area  to  bay  area. 

The  idealized  estuary  shape  is  approximated  by  a  semi-ellipse  with  two  500-m-wide  sand  bars 
and  a  single  tidal  inlet  between  them  (Figure  2).  The  estuary  is  medium  sized  with  a  total  area  of 
40  km“.  The  offshore  bathymetry  follows  an  equilibrium  profile  for  a  median  sediment  size  of 
0.1  mm.  The  inlet  length  and  width  are  500  and  220  m,  respectively.  The  tidal  amplitude  is  0.5  m 
(micro-tidal)  and  the  bay  depth  is  4  m  (relative  tidal  amplitude  =  0.125). 

Seelig  and  Sorensen  (1978),  Boon  (1975),  and  Boon  and  Byme  (1981)  showed  that  the  rate  at 
which  the  wetted  basin  surface  area  varies  with  depth  or  hypsometric  variation  is  a  key  factor 
controlling  tidal  flows  and  tidal  asymmetry.  The  wetted  surface  area  as  a  function  of  the  water 
depth  is  given  by  the  idealized  curve  (Boon  1975), 


r_  r(l-G) 
r  +  G(l  —  r) 


(1) 


in  which  A'  =  (A  -  Amin)/(Amax  -  Amm),  A  is  the  wetted  surface  area  at  water  elevation  h,  and  Amax  is 
the  wetted  surface  area  at  water  depth  hmax.  The  fitted  coefficient  r  =  AmJAmax,  where  Amin  is  the 
minimum  area  at  relative  height.  The  parameter  G  =  (1  -  h')y,  in  which  y  is  a  positive  exponent 
(shape  factor)  that  determines  the  shape  of  the  curve,  and  h'  =  (h  -  hmin)/(hmi_lx  -  /?mm),  where  hmm 
and  hmax  are  the  minimum  and  maximum  values  of  h  (i.e.,  /?mm  =  mean  low  water  (MLW)  and 
hmax  =  mean  high  water  (MHW)).  Figure  3  shows  how  the  Equation  1  varies  as  a  function  of  r 
and  y. 
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Figure  2.  Idealized  backbarrier  estuary  (r  =  0. 1 5  and  y  =  2).  Depth  is  relative  to  mean  sea 
level,  and  square  cells  represent  inactive  regions. 


Figure  3.  Idealized  estuary  area-height  curves. 
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Below  the  intertidal  range,  the  estuary  hypsometric  curve  gradually  approaches  zero.  The  values 
for  Am in  and  r  were  estimated  as  11.8  km2  and  0.15,  respectively,  by  comparing  the  percent 
wetland  coverage  and  ratio  of  intertidal  to  estuarine  channel  water  volume.  These  estimates  are 
consistent  with  the  values  reported  in  Friedrichs  and  Aubrey  (1988)  and  references  therein.  For 
r  =  0.15,  the  shape  of  the  hypsometric  curve  is  approximately  convex  for  y  <  2  and  becomes 
concave  for  y  >  2 . 5 . 

The  change  in  hypsometric  curve  due  to  wetland  loss  is  simulated  by  increasing  the  shape  factor 
y.  Five  numerical  tests  were  conducted  with  shape  factors  equal  to  1,  2,  4,  5,  and  6.  Because 
vegetation  can  grow  only  within  a  portion  of  the  intertidal  range,  the  amount  of  habitat  area  for 
vegetation  decreases  as  the  shape  factor  increases.  As  an  example,  Figure  3  shows  the 
bathymetry  for  an  idealized  estuary  with  y  =  2.  The  fundamental  assumptions  here  are  (1)  the 
inlet,  estuary  bottom,  and  perimeter  do  not  change  (r  =  constant),  and  (2)  for  varying  degrees  of 
wetland  loss,  the  area-depth  relationship  of  the  estuary  within  the  intertidal  zone  can  be 
approximated  by  Equation  1 . 

Stable  intertidal  salt  marshes  occupy  relatively  flat  areas  known  as  marsh  platforms  within  the 
intertidal  zone.  The  marsh  vegetation  modeled  herein  is  intended  to  represent  Spartina 
altemiflora  or  smooth  cord  grass,  which  is  the  most  prevalent  salt  marsh  vegetation  in  the  United 
States.  The  range  of  S.  altemiflora  depends  mostly  on  the  mean  range  of  tide  (Mn)  defined  as 
Mn  =  MHHW  -  MLLW.  For  the  modeling  discussed  herein,  the  upper  and  lower  elevation  limits 
of  vegetation  are  calculated  as  zupper  =  0.43Mn  +  0.25,  and  ziower  =  -0.24Mn  +  0.09,  respectively, 
from  Mckee  and  Patrick  (1988).  As  an  example,  the  vegetation  distribution  for  the  idealized 
estuaries  with  y  =  2  is  shown  in  Figure  4.  Stem  densities  of  S.  altemiflora  vary  both  spatially 
within  a  marsh  and  from  region  to  region  and  range  from  16  to  4,080  stems/m2,  but  are  typically 
in  the  range  of  50  to  350  stems/m"  (Visser  et  al.  1996  and  references  therein).  Here,  a  relatively 
dense  salt  marsh  density  of  320  stems/m2  is  used  with  an  average  stem  diameter  of  0.2  cm  and  is 
constant  over  the  vegetated  salt  marsh  (Leonard  and  Luther  1995;  Lightbody  and  Nepf  2006). 

Vegetation  Flow  Drag.  The  estimation  of  vegetation  flow  drag  is  not  straightforward  in  a 
regional  model  because  of  the  variation  in  plant  density  and  micro-topography.  First-  and 
second-order  tidal  channels  within  the  marsh  have  the  effect  of  reducing  flow  drag.  The 
vegetation  flow  drag  usually  requires  field  measurements  for  calibration.  Vegetation  flow  drag 
should  be  represented  in  a  physically  reasonable  form  and  not  by  simply  increasing  the  bottom 
friction.  Bottom  friction  decreases  as  the  water  depth  increases,  whereas  the  opposite  is  true  for 
vegetation  flow  drag.  Here,  vegetation  is  modeled  as  rigid  vertical  cylinders,  and  a  correction  is 
applied  to  compensate  for  the  presence  of  tidal  channels.  Vegetation  flow  drag  is  estimated  as 
(Wu  and  Wang  2004) 


x  =  a  —  C\p  X  lulu 

v  c  2  dr  w  a  \  | 


(2) 
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Figure  4.  Example  vegetation  coverage  for  an  idealized  backbarrier  estuary  (r  =  0.15 
and  y  =  2).  Blue  areas  represent  a  vegetation  density  of  320  stems/m2  and  square  cells 

represent  inactive  regions. 


where  ac  is  a  correction  factor  for  the  presence  of  a  tidal  channel  network  (approximated  here  as 
the  ratio  of  vegetated  marsh  platform  area  to  total  marsh  area),  Cd  is  an  empirical  bulk  drag 
coefficient  (set  to  1.0),  pw  is  the  water  density,  Xa  is  the  vegetation  projection  area  per  unit 
volume,  and  U  is  the  flow  velocity  vector.  The  projected  vegetation  area  per  unit  volume  is 
calculated  as  Xa  =  nvDv  mini/?,.,  h)  where  nv  is  the  number  of  stems  per  square  meter,  Dv  is  the 
stem  diameter,  hv  is  the  average  plant  height  (set  to  0.75  m),  and  h  is  the  local  water  depth. 
Equation  2  enters  the  momentum  equations  in  a  similar  manner  as  the  bottom  friction. 

Model  Setup.  The  idealized  estuary  was  forced  with  the  principal  lunar  semidiurnal  constituent 
M2  with  an  amplitude  of  0.5  m.  For  a  tidal  range  of  1  m,  the  lower  and  upper  limits  of  vegetation 
are  -0.15  m  and  0.68  m  with  respect  to  mean  sea  level.  The  grid  resolution  varied  from  25  m  at 
the  tidal  inlet  to  50  m  elsewhere.  A  constant  Manning’s  n  of  0.025  was  specified  over  the 
calculation  domain. 

Numerical  Model  Results.  Example  current  velocities  at  four  observation  stations  are  shown 
in  Figure  5  for  the  idealized  estuaries  with  shape  factors  equal  to  2  and  5,  which  represent  vege¬ 
tated  areas  of  69  and  36  percent  (48  percent  wetland  loss).  The  results  show  an  increase  in  both 
maximum  flood  and  ebb  current  velocities  at  the  tidal  inlet  of  25  and  21  percent,  respectively. 
Increased  current  velocities  are  also  observed  at  all  other  stations  with  ranges  between  30  and 
46  percent.  The  change  in  sediment  transport  rate  would  be  significantly  greater  than  the  current 
velocity. 
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Figure  5.  Current  velocities  in  the  x  (flood)  direction  at  four  stations  for  the  idealized  inlet: 

(a)  shape  factor  =  2,  and  (b)  shape  factor  =  5. 


Key  numerical  model  parameters  and  a  summary  of  results  are  listed  in  Table  1.  The  proportion 
of  the  bay  that  is  covered  by  vegetation  varies  between  78  and  26  percent.  The  duration  of  the 
flood  and  ebb  tides  and  the  ratio  of  the  maximum  flood  and  ebb  currents  are  relative  indicators  of 
flood  and  ebb  dominance.  The  increase  in  tidal  prism  with  increasing  wetland  loss  causes  an 
increase  in  both  the  flood  and  ebb  current  velocities  through  the  tidal  inlet.  However,  as  the 
shape  factor  increases,  the  bathymetric  profile  becomes  more  concave,  reducing  the  ebb  domi¬ 
nance.  Although  all  of  the  cases  tested  were  ebb  dominant,  represented  by  the  ratio  of  the  maxi¬ 
mum  flood  to  ebb  current  velocity,  there  is  a  clear  tendency  toward  flood  dominance  with 
increasing  shape  factor.  The  cases  show  only  a  slight  increase  (decrease)  in  flood  (ebb)  tides  of 
0.25  hr. 
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Table  1 

Description  of  five  idealized  estuaries  with  shape  factors  between  1  and  6  and 
summary  of  current  velocities  at  the  tidal  inlet  (nvDv  =  0.64  stem/m). 


Parameter 

Case  A 

Case  B 

Case  C 

Case  D 

Case  E  j 

Shape  factor  y 

1 

2 

4 

5 

6 

100  x  AJAb,  % 

78 

69 

47 

36 

26 

Tidal  prism,  km3 

48.03 

60.02 

72.73 

76.17 

78.69 

V\Nb 

0.79 

1.17 

1.48 

1.52 

1.52 

Flood  tide,  hr 

6.75 

6.75 

6.75 

7.00 

7.00 

Ebb  tide,  hr 

5.25 

5.25 

5.25 

5.00 

5.00 

^tiood  >  m/sec 

0.86 

1.10 

1.38 

1.45 

1.50 

Ue bb,  m/sec 

1.10 

1.33 

1.60 

1.68 

1.73 

l^flood/^ebbl 

0.79 

0.82 

0.86 

0.86 

0.87  ; 

NOTE: 

V  =  Wetland  intertidal  water  volume  of  wetland. 
Vb  =  Estuarine  water  volume. 

Ufiood  =  Maximum  flood  current  velocity  at  inlet. 
Ue bb  =  Maximum  ebb  current  velocity  at  inlet. 

A,  =  Wetland  intertidal  surface  area. 

Ab  =  Estuarine  surface  area. 


DISCUSSION:  Recognizing  possible  limits  of  the  model  assumptions,  including  the  idealized 
distribution  of  estuarine  vegetation,  hypsometric  curve,  and  simplified  tidal  forcing,  the  results  of 
the  idealized  estuary  simulations  provide  quantitative  insight  into  how  wetland  loss  and  restor¬ 
ation  in  backbarrier  bays  modify  current  velocities  at  tidal  inlets.  Two  main  processes  are 
observed.  First,  as  the  tidal  prism  increases  due  to  wetland  loss,  there  is  an  increase  in  both  the 
flood  and  ebb  current  velocities.  Second,  as  wetland  loss  occurs  and  the  bay  increases  in  depth, 
the  bathymetric  profile  becomes  more  convex,  decreasing  ebb  dominance.  As  the  estuary 
becomes  less  ebb  dominant,  less  sediment  is  exported,  thereby  preventing  the  estuary  from 
eroding.  In  a  similar  manner,  as  wetlands  are  restored  or  developed  naturally,  the  intertidal  zone 
increases  and  ebb  dominance  is  enhanced,  thus  increasing  the  export  of  sediment  from  the 
estuary.  This  interaction  between  estuary  geometry  and  tidal  asymmetry  is  a  natural  way  by 
which  estuaries  stabilize  themselves. 

Actual  estuaries  have  more  complex  configurations  and  forcing  than  the  idealized  estuaries 
presented  here.  Other  factors  are  involved  in  the  morphologic  evolution  of  estuaries  such  as 
sediment  inputs  and  relative  sea  level  rise.  For  estuaries  to  be  stable  in  a  long-term  sense  they 
must  keep  up  with  relative  sea  level  rise  by  accreting  both  mineral  and  organic  material.  In  the 
case  of  Barataria  Bay,  eustatic  sea  level  rise,  erosion,  and  local  subsidence  are  drowning  the 
wetlands  in  place.  As  bay  area  increases,  fetch  increases  and  wind-generated  waves  cause  marsh 
fringe  erosion. 

In  this  study,  the  vegetation  density  was  assumed  constant  with  respect  to  tidal  elevation,  but  in 
nature  this  is  not  the  case.  Vegetation  occurs  most  commonly  at  specific  elevations  (marsh 
platform).  The  predominant  marsh  elevation  will  vary  for  each  location,  depending  on  the  local 
physical,  chemical,  and  weather  conditions.  It  is  expected  that  greater  vegetation  densities  in 
shallow  water  will  enhance  ebb  dominance. 
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The  idealized  modeling  approach  described  here  did  not  include  the  interactions  between 
morphology  change  and  hydrodynamics  such  as  shore  erosion,  changes  in  tidal  inlet  cross- 
sectional  area,  and  evolution  of  ebb  and  flood  shoals.  Only  the  resulting  change  in 
hydrodynamics  was  considered  for  different  conditions  representing  various  degrees  of  wetland 
coverage.  The  inlet  and  bay  geometry  were  held  constant,  while  the  change  in  shape  of  the 
intertidal  zone  for  varying  degrees  of  wetland  loss  was  simulated  using  an  idealized  curve.  It  is 
unclear  how  different  causes  of  wetland  loss  such  as  drowning  and  marsh  fringe  erosion  produce 
different  responses  in  the  hypsometric  curve.  Additional  simulations  are  needed  for  various  bay 
sizes,  inlet  dimensions,  and  tidal  regimes  to  obtain  a  more  complete  understanding  of  the 
interactions  between  wetlands  and  tidal  inlets. 

The  following  areas  were  identified  for  further  study  in  model  development: 

1 .  Vertical  distribution  of  wetland  vegetation  and  how  it  changes  tidal  asymmetry. 

2.  Change  in  estuary  bathymetry  and  geometry  (area-height  relationship)  caused  by  wetland 
loss. 

3.  Interactions  between  long-term  morphology  change  and  hydrodynamics. 

4.  Influence  of  tidal  regime  and  relative  tidal  amplitude  to  estuary  resilience  to  relative  sea 
level  rise  and  wetland  loss. 

5.  Influence  of  stonns  and  atmospheric  fronts  on  the  overall  morphological  evolution  of 
estuaries  (will  vary  regionally). 

Because  flow  resistance  caused  by  vegetation  is  proportional  to  the  velocity  squared,  wetland 
vegetation  is  most  efficient  at  dissipating  momentum  in  areas  of  stronger  currents.  This  leads  to 
considerations  about  the  location  of  vegetation  and  mud  flats  within  estuaries  and  wetlands. 
Consequently,  it  is  hypothesized  that  constructing  and  restoring  vegetated  marshes  in  areas  of 
relatively  weak  current  velocities  and  mud  flats  in  regions  of  relatively  strong  currents  (usually 
occurring  near  tidal  inlets)  should  promote  flood  dominance.  These  placement  alternatives  must 
be  carefully  evaluated  to  ensure  that  the  hydraulic  characteristics  are  adequate  for  the 
sustainability  of  the  mud  flats  and  wetlands,  and  waves  and  currents  will  not  erode  the  project 
site.  Tidal  theory  suggests  that  wetland  channels  that  are  more  shallow  and  narrow  while 
maintaining  the  same  cross-sectional  area  will  promote  flood  dominance  through  non-linear 
frictional  effects.  It  is  recommended  that  the  functional  design  of  constructed  wetlands  be 
supported  by  hydrodynamic,  sediment  transport,  and  particle-tracking  models. 

CONCLUSIONS:  The  long-term  progression  of  tidal  inlets  and  navigation  channels  is  closely 
related  to  the  evolution  of  estuarine  wetlands  through  the  tidal  prism,  sediment  exchange 
between  the  open  ocean  and  estuary,  the  estuary  surface  area  as  a  function  of  bay  area  and  depth, 
vegetative  damping,  and  sediment  balance.  Wetland  loss  is  expected  to  increase  in  many 
estuaries  due  to  sea  level  rise,  climate  change,  and  reduced  sediment  supply.  The  consequences 
of  wetland  loss  on  the  tidal  inlet  morphology  have  been  observed  over  the  past  100  years  in 
Barataria  Bay,  LA.  A  numerical  experiment  of  an  idealized  estuary  demonstrates  the  increase  in 
both  flood  and  ebb  current  velocities  and  the  tendency  to  shift  toward  flood  dominance  with 
increasing  wetland  loss.  Future  work  will  provide  predictive  tools  and  guidance  on  the  design 
and  maintenance  of  tidal  inlets  and  navigation  channels  that  can  take  into  account  the 
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interactions  among  wetland,  tidal  inlets,  and  navigation  channels.  New  tools  for  studying 
estuarine-wetland  interactions  are  under  development,  and  user  input  is  requested  on  modeling 
capability  needs. 

ADDITIONAL  INFORMATION:  This  CHETN  is  a  product  of  the  Inlet  Channels  and  Estuaries 
Work  Unit  of  the  Coastal  Inlets  Research  Program  being  conducted  at  the  U.S.  Army  Engineer 
Research  and  Development  Center,  Coastal  and  Hydraulics  Laboratory,  Vicksburg,  MS. 
Questions  about  this  technical  note  can  be  addressed  to  Alejandro  Sanchez  (Voice:  601-634- 
2027,  Fax:  601-634-3433,  e-mail:  Alejandro. Sanchez@usace.army.min.  For  information  about 
the  CIRP,  please  contact  the  CIRP  Program  Manager,  Dr.  Nicholas  C.  Kraus  (Voice:  601-634- 
2016,  e-mail:  Nicholas.C.Kraus@usace.army.mil).  The  CHETN  benefited  from  technical  reviews 
by  Julie  Dean  Rosati,  Dr.  Kraus,  and  Sophie  Munger.  This  document  should  be  cited  as  follows: 

Sanchez,  A.  2008.  Interactions  between  wetlands  and  tidal  inlets.  Coastal  and 
Hydraulics  Engineering  Technical  Note.  ERDC/CHL  CHETN-IV-72.  Vicksburg, 

MS:  U.S.  Army  Engineer  Research  and  Development  Center. 
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